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Reversible diffusion limited cluster aggregation of hard spheres with rigid bonds was simulated 
and the self diffusion coefficient was determined for equilibrated systems. The effect of increasing 
attraction strength was determined for systems at different volume fractions and different interaction 
ranges. It was found that the slowing down of the diffusion coefficient due to crowding is decoupled 
from that due to cluster formation. The diffusion coefficient could be calculated from the cluster size 
distribution and became zero only at infinite attraction strength when permanent gels are formed. 
It is concluded that so-called attractive glasses are not formed at finite interaction strength. 

PACS numbers: 05.10.Ln, 82.70.Dd, 82.70.Gg 



I. INTRODUCTION 

Reversible aggregation of small particles in solution 
is a common phenomenon. It leads to different equi- 
librium states depending on the volume fraction (0) of 
the particles and the strength of the interaction energy 
(m). Weak attraction results in the formation of tran- 
sient aggregates at low </> and a transient percolating 
network at high 0, while strong attraction may drive 
phase separation into a high and a low density liquid 



10 



The strength of the interaction and thus the equilib- 
rium properties are determined by the ratio of the bond 
formation ( a) and the bond brea king {(3 ) probab ility 

P, H [si, IM [13, ii, IM EO, IM El m S, El , 

while the kinetics of such systems depend on the abso- 
lute values of a and (3. Two limiting cases may be dis- 
tinguished: diffusion limited cluster aggregation (DLCA) 
for which a bond is formed at each collision (a = 1) and 
reaction limited cluster aggregation (RLCA) for which 
the probability to form bonds goes to zero (a ^ 0). 

The average long time self diffusion coefficient {Di) 
of non-interacting hard spheres decreases with increas- 
ing volume fraction [i^, [4^, [1^ [1^, [s^l and the system 
forms a glass above a volume fraction of about 0.585. 
For a recent review of theories and experiments on the 
glass transition in colloids see Sciortino and TartagliafsTll. 
Mode coupling theory predicts that Di goes to zero at a 
critical volume fraction 0c = 0.516 following a power law: 
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The pictorial view is that particles become trapped in 
cages formed by neighbourin g! p articles. Experiments 
[5^ and computer simulations jl9l. [ssj confirmed this be- 
haviour over a range of (j) though with a larger value 
of (f)c (0.585). However, the diffusion coefficient is not 
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truly zero at (pc and some mobility is still possible be- 
tween 0c and the volume fraction of random close packing 
{'t'cp = 0.64). The reason is that fluctuations of the cage 
size around the average value, allow particles to "hop" 
from one cage to another. Mode coupling theory uses 
the static structure factor as input and therefore cannot 
include the effect of fluctuations and heterogeneity. 

Introducing reversible bond formation between the 
hard spheres has two consequences. Firstly, the struc- 
ture factor is modified, because on average more parti- 
cles will be within each others bond range. Secondly, the 
diffusion of bound particles becomes correlated for some 
duration that is related to the bond life-time. Obviously, 
the latter effect leads to a decrease of Di, because clus- 
ters of bound particles move more slowly than individual 
spheres and not at all in the case when they form a per- 
colating network. However, modification of the structure 
may also lead to an increase of Di compared to the hard 
sphere case [H, [l^, [13], because the accessible volume 
for a particle in which it can move increases. In other 
words, the average cage size increases. This effect can be 
clearly seen for the self diffusion of tracer particles in a 
medium of fixed spherical obstacles. For a given volume 
fraction of obstacles, Di is smaller when the obstacles are 
randomly distributed than when they form a percolating 
network [s^ . 

An increase of Di with increasing attraction has actu- 
ally been observed for concentrated suspensions of hard 
spheres with weak short range attraction where the par- 
ticles can freely move within the attraction range. With 
further increase of the interaction energy the effect on 
the structure weakens, but the bond life-time increases 
so that Di decreases again. These effects have also 
been found in molecular dynamics simulations of such 
systems [H, [l^ [sj- Mode coupling theory describes 
this behaviour purely on the basis of the changes in the 
static structure factor and predicts complete arrest at 
finite interaction strength even at low volume fractions 
[11]. Systems arrested by attraction are called attrac- 
tive glasses [13] to distinguish them from the ordinary 
repulsive glasses f52j . 

Obviously, irreversible aggregation, i.e. infinite attrac- 
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tion strength, leads to arrest at any volume fraction when 
all spheres become part of the percolating network [sl] . 
Generally, arrested systems at low volume fractions due 
to irreversible aggregation are called gels. In analogy, 
percolated systems with finite interaction strength may 
be called transient gels, though in practice one would only 
do that only if the bond life-time is very long so that the 
system flows very slowly. The issue addressed here is how 
Di varies with increasing attraction strength and specif- 
ically whether Di becomes zero in transient gels with a 
finite bond life-time. In other words, whether transient 
gels can be attractive glasses. Even if, as for repulsive 
glasses, Di does not truly become zero due to hopping 
processes, the question remains whether Di can be use- 
fully described as going to zero at a finite critical inter- 
action energy following a power law over a broad range 
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as was suggested in the literature (54|, |59| (uc). 

As mentioned above, phase separation occurs when the 
attraction is strong, and needs to be avoided if one wishes 
to study the effect of strong attraction on Di at low vol- 
ume fractions. This can be done to some extent by lim- 
itin g th e maximum number of bound neighbours below 6 
[s^ lSOll or by introducing a correlation between the bond 
angles with different neighbours [6]| . One can also in- 
troduce a long range repulsive interaction to push phase 
separation to lower concentrations [13, [H, [6^. Phase 
separation disappears completely if the bonds are rigid 
and the bond range is zero, because only binary colli- 
sions can occur so that the average coordination number 
cannot exceed two [l^ . For this model Di can be studied 
at any interaction strength and volume fraction without 
the interference of phase separation. 

An extensive study of the structural properties of such 
systems was published earlier [11] . It was shown that the 
equilibrium properties depend on the ratio a/ (3 and the 
elementary step size s of the Brownian motion. As long 
as s is much smaller than the average distance between 
the particles, the equilibrium properties are determined 
by a single parameter called the escape time (tg) that is 
a combination of s, a and /?. is defined as the ex- 
cess time it takes for two particles at contact to become 
decorrelated compared to the situation where no bonds 
are formed (a = 0). If the interaction range is finite and 
if s is much smaller than the interaction range, then the 
escape time becomes independent of s and the equilib- 
rium properties are the same as for spheres interacting 
with a square- well potential. 

Here we present results of the mean square displace- 
ment (MSD) of the particles as a function of the interac- 
tion strength for systems forming rigid bonds with zero 
and finite interaction strength. The main conclusion of 
this work is that introducing attraction leads to a de- 
crease of Di without any sign of arrest at finite interac- 
tion strength. As mentioned above, there are two differ- 
ent mechanisms for slowing down: one is controlled by 
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FIG. 1: MSD of spheres starting from a random position at 
cj> = 0.49 for dilTerent te as indicated in the figure. 



the bond life-time and the other is controlled by crowd- 
ing. We show here that for rigid bond formation with 
a — 1 the two mechanisms are independent and can be 
factorized. Since the first mechanism leads to complete 
arrest only for irreversible aggregation, we propose that 
the expression "attractive glass" is abandoned in favour 
of the expression "gel" that is commonly used to describe 
irreversibly bound percolating networks. The results pre- 
sented here resemble to some extent those obtained re- 
cently by Zaccarelli et al.[59l| for attractive spheres with 
limited valency. The main difference is that in those 
simulations the bonds were not rigid and localized mo- 
tion was possible even at infinite bond strength. We will 
briefly discuss the effect of bond flexibility. 



II. SIMULATION METHOD 

The simulation method for aggregation with zero in- 
teraction range has already been detailed elsewhere [i^ . 
Initially, TV non-overlapping spheres with unit diameter 
are positioned randomly in a box with size L, using peri- 
odic boundary conditions, so that — {tt/6)N/L'^. Dif- 
ferent sizes between L = 10 and L = 50 were used, and 
the results shown here are not influenced significantly by 
finite size effects. N times a particle is chosen and moved 
a distance s in a random direction. When the movement 
leads to overlap with another sphere it is truncated at 
contact. After this movement step all spheres in contact 
are bound with probability a leading to the formation of 
Nc clusters that are defined as sets of bound particles. 
In the next movement step, Nc times a cluster is cho- 
sen and moved a distance s in a random direction with 
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FIG. 2: Dependence of the self diffusion coefficient of ran- 
domly distributed hard spheres on the volume fraction ob- 
tained from computer simulations: present work (circles), [1^ 
(solid line), |48|] (squares) and [49(| (triangles). 

a probability that is inversely proportional to the diam- 
eter of the cluster, thus simulating non-draining within 
the clusters. In the following cluster construction step, 
all bound particles are broken with probability (3 and 
bonds are formed for new contacts with probability a. 
Movement and cluster formation steps are repeated until 
all particles have formed a single cluster for irreversible 
aggregation or until equilibrium is reached for reversible 
aggregation. 

The diffusion coefficient was calculated as (r^)/(6i) 
and the time unit was defined as the time needed for a 
single particle to diffuse its own diameter so that the dif- 
fusion coefficient at infinite dilution is Z)o = 1/6- In this 
article we will present diffusion coefficients normalized by 
Dq. For the case of irreversible (/3 — 0) diffusion limited 
(a = 1) aggregation, the kinetics of cluster growth were 
the same as predicted from the Smolechowski equations 
for DLCA if s was chosen sufficiently smallfs^. As men- 
tioned in the introduction, the equilibrium properties for 
reversible aggregation with zero interaction range are de- 
termined by the escape time psj : 

<e = (2.37s-|-s2)^ (3) 

At equilibrium, a distribution of self-similar transient 
clusters is formed together with a transient percolating 
network if and 4> are sufficiently large. 

Reversible aggregation with finite interaction range e 
was simulated by forming bonds with probability a when 
the centre to centre distance between particles is less than 
1 -I- e [1^. Bonds were again broken with probability 
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FIG. 3: MSD of spheres during irreversible aggregation 
(DLCA) for different volume fractions as indicated in the fig- 
ure. 

/?. In this case the equilibrium properties are indepen- 
dent of the step size if s e. The system is equivalent 
to particles interacting through a square well attraction 
with interaction energy u = ln(l — P), where P is the 
probability that particles within each others interaction 
range are bound: P — a/(a. + (3). m is given in units 
of the thermal energy and is equivalent to the inverse 
temperature that is sometimes used to express the inter- 
action strength. The equilibrium properties of systems 
with different interaction range are close if they are com- 
pared at the same second virial coefficient (52) [M]- ^2 is 
the sum of the excluded volume repulsion and the square 
well attraction: B2 = B^ep — Batt- Brep = 4 and Batt 
is determined both by the interaction strength and the 
interaction range fssj : 

Batt - 4[(1 + e)3 - 1] [exp (-u) - 1] (4) 

in units of the particle volume 



III. RESULTS 

A. Zero interaction range 

Fig. [1] shows the mean square displacement, i.e. (r^), 
of spheres at (j) — 0.49 for different values of <e with 
a — 1, starting from a random distribution. Initially, the 
particles move freely until they begin colliding with other 
particles. This leads to slowing down of the MSD even in 
the absence of attraction {te =0). At long times (r^) is 
again proportional to t, but the long time diffusion coef- 
ficient (Di) is smaller. Fig. [5] shows that the dependence 
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FIG. 4: MSD for equillibriated systems (a) at (/> = 0.49 for 
different te. Fig. 4b shows a mastercurve obtained by plotting 
the same data as a function of Dst. 



of the normalized diffusion coefRcient of hard spheres in 
the absence of attraction (Dhs) on is in good agreement 
with Hterature results O, [13, [H, O, 50]. 

Introducing attraction slows down the MSD, because 
transient clusters are formed that move more slowly and 
above a critical value of a transient percolating net- 
work is formed. Particles that are part of the transient 
network are immobile and need to break bonds before 
they can diffuse. In the limit of = oo, correspond- 
ing to irreversible DLCA, all particles are permanently 
stuck when they become part of the percolating network. 
Naturally, the value of (r^) where the particles get stuck 
decreases with increasing volume fraction, see Fig [31 

Once the system has reached equilibrium one can start 



again measuring the MSD. Fig. SJj shows the MSD for 
the same systems as in Fig[Tl but starting from the equi- 
librium state. In this case the initial diffusion is slower 
than that of individual spheres, since the system contains 
clusters and for te > 0.43 a percolating network. After 
some time the clusters collide causing a further slowing 
down of the MSD, but simultaneously, for a = 1, bonds 
are formed and broken. Of course, for t ^ the same 
long time diffusion coefScient is obtained independent of 
the starting configuration, compare Fig. [1] Remarkably, 
the effect of particle collisions on the slowing down of 
the MSD is the same at different te- Therefore the MSD 
obtained at different can be superimposed within the 
statistical error by simple time shifts (see Fig. Ht). The 
implication is that Di can be factorized in terms of the 
short time diffusion coefficient {Ds) and Dhs- 
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At short times, i.e. before bond breaking becomes sig- 
nificant and before the particles collide, the MSD is de- 
termined by free diffusion of the clusters. Consequently, 
Ds can be calculated from the size distribution of the 
clusters {N{m)) which is a function of and 
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where D{m) is the average free diffusion coefScient of 
clusters with aggregation number to. If a percolating 
network is present we need to consider only the sol frac- 
tion of mobile clusters (F): 



F_ TO ■ N{m) ■ D{m) 
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where m„ is the number average aggregation number of 
the sol fraction. In the absence of hydrodynamic inter- 
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For the more 



realistic case of non-draining clusters Dim) is inversely 
proportional to the radius of the clusters so that Dg de- 
creases more weakly with increasing to„. 

Figure [H] shows Dg as a function of tf. for (f) = 0.49 and 
(f) = 0.05. We have calculated Dg from the cluster size 
distribution using Eq. ([6]) and Di was calculated using 
Eq. ([5]). Comparison with the simulation results, see 
Fig I®, demonstrates that the effect of attraction on Dg 
is indeed fully determined by the cluster size distribu- 
tion. Initially, Di decreases with increasing ig, because 
TO„ increases until at a critical value of (^e) ^ Per- 
colating network is formed of immobile particles, i* is 
0.43 and 180 for (j) = 0.49 and (f) = 0.05, repectively. At 
the percolation threshold the weight average aggregation 
number diverges (45j . but to„ and thus Di remain finite. 
The maximum value of to„ is obtained at the percolation 
threshold, beyond which it decreases with increasing vol- 
ume fraction. Beyond the percolation threshold the sol 
fraction decreases much more strongly with increasing at- 
traction than TO„ so that Di continues to decrease. How- 
ever, it is obvious that Di will become zero only if = 0, 
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FIG. 5: The short time diffusion coefficient as a fuction of te 
for (f) = 0.49 (squares) and 4> = 0.05 (circles) . Closed and 
open symbols represent simulation results and calculations 
using Eq. ([6}. 
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FIG. 7: Probability distribution of at different times 0.001 
(circle), 0.3 (triangledown), 0.5 (square), 3 (diamond) and 167 
(triangle) for (f} — 0.49 and — 2. The solid lines represent 
the distributions for freely diffusing spheres with the same 
average MSD, see Eq. |8l 



1.0 



0.8 - J 
0.6 - / 

s - / 
s : J 

0.4 - J 
0.2 - / 

0.0 ^ 

10-5 10-* 10-3 10-2 10-1 IQO 

t 

FIG. 6: Fraction of mobile particles as a fuction of time for 
an equillibrated system at te = 2 and — 0.49 . The solid 
line is a guide to the eye. 



i.e. at infinite attraction strength. At large values of t^, 
Di decreases linearly with increasing ig. As expected, at 
lower volume fractions stronger attraction is needed to 
cause significant slowing down. 

The MSD shown in Fig. 0] represent an average over 
all particles, and do not show how the displacement of 



the particles is distributed. In the absence of attraction 
the probability distribution that a particle has moved a 
distance at a given time (P(r^)) is given by: 

^('^''^^-(i^^^^P^-i^) 

because each particle is equivalent. However, in the pres- 
ence of attraction the displacement is highly heteroge- 
neous, because the particles belong to clusters of differ- 
ent sizes. We have illustrated this for (j) = 0.49 at = 2. 
In this case almost all of the particles belong to the gel 
fraction and cannot move until their bonds are broken. 
Thus at short times the displacement is highly heteroge- 
neous with a large fraction of particles that do not move 
at all and a small fraction of sol particles (F) that dif- 
fuse freely until they collide. With increasing time, more 
and more of the sol particles collide with the percolating 
network and stop moving, while more and more gel par- 
ticles break lose and start diffusing. Fig [S] shows how the 
fraction of mobile particles that has moved after a time 
t {M{t)) increases with increasing t. 

Figure [7] compares P{r'^) at different times with the 
distribution that would have been found for same (r^) 
if all particles had been equivalent. Note that the area 
under the curves is equal to M{t). The heterogeneity of 
the displacement decreases with increasing time until for 
t ^ te all particles have formed and broken bonds many 
times so that M{t) = 1 and -P(r^) is given by Eq. [8] . 
Similar observations were made at 4> = 0.05. 

Figure [7] resembles the results presented by Puertas et 
al. [66j . However, they used molecular dynamics simu- 
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FIG. 8: The short time diffusion coefficient at = 0.49 (cir- 
cles), (j> ~ 0.3 (squares) and = 0.05 (triangles) for two dif- 
ferent interaction ranges e = 0.5 (closed) and e = 0.1 (open) 
as a function of Ban- The solid lines are guide to the eye. 

lations for which bound particles move freely as long as 
the displacement does not involve breaking of bonds, as 
will be discussed below. For this reason they observed a 
peak situated at a small r that represents the displace- 
ment of bound particles. A second peak was found at 
larger values of r representing the sol fraction and moved 
to longer distances with increasing time. The amplitude 
of this peak increased with increasing time as more and 
more particles break their bonds. 

B. Finite interaction range 

The effect of finite interaction range was tested for 
e = 0.1 and e = 0.5. We found also for this case that 
Di could be factorized into Dg and Dhs, see Eq. [SI 
and that the effect of attraction on Di is fully deter- 
mined by the free diffusion of the clusters. However, the 
dynamic heterogeneity that was important for the zero 
range interaction disappeared rapidly with increasing in- 
teraction range especially at higher concentrations, be- 
cause monomers and clusters can move only a short dis- 
tance before they interact. Individual particles exchange 
rapidly between different clusters including the gel frac- 
tion, at least for a = 1. Therefore each particle explores 
more rapidly the different dynamics of the system. 

Recently, it was shown that the cluster distribution in 
equilibrium is similar for the two interaction ranges if 
compared at the same values of the second virial coeffi- 
cient, at least for (j) up to 0.2 [2^. We therefore expect 
that Dg is the same at the same i?2- Fig. [5] shows Dg 
as a function of Batt for 3 different volume fractions with 
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FIG. 9: The short diffusion coefficient as a function m„/_F for 
different volume fractions {(j)) and interaction ranges (e) as 
indicated in the figure . The solid line represents Ds — mn/F 

e — 0.5 and e = 0.1. As mentioned in the introduction, 
for lower volume fractions the slowing down of Dg can be 
studied only over a limited range before phase separation 
occurs. As expected, Dg decreases with increasing Batt 
and the decrease is more important at higher concentra- 
tions. Again no sign of critical slowing down of the MSD 
was observed. For e = 0.5 and e = 0.1 the dependence 
of Dg on Batt was similar for (f) — 0.05 and 0.3, but for 
(f) = 0.49 it was stronger for e = 0.1 than for e — 0.5. At 
high volume fraction B2 is no longer the main parame- 
ter that determines the cluster size distribution because 
higher order interaction becomes important. 

Recently, Foffi. et al. [1^ reported a simulation study of 
the effect of the interaction range on Di for hard spheres 
with a square well interaction. Molecular dynamics sim- 
ulations were used that gave the same equilibrium struc- 
tures as with the method used here. However, molecu- 
lar dynamics simulations gives different dynamics so that 
absolute values of Di cannot be compared with the re- 
sults presented here. Nevertheless, they also found that 
the MSD of particles was independent of the interaction 
range if compared at the same value of B2 and they ar- 
gued that Di was the same because the number of bonds 
and the bond life-time was determined by B2. We expect 
that deviations will be found also with this method at 
higher concentrations and larger interaction ranges when 
higher order interaction becomes important. 

IV. DISCUSSION 

There are two causes for the decrease of Di in a system 
of attractive hard spheres. The first one is collisions with 
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other spheres. This effect of crowding increases with in- 
creasing volume fraction and leads to strong decrease of 
Di near (p « 0.58 that can be described to some extent 
by mode coupling theory, as mentioned in the Introduc- 
tion. The second cause is bond formation, which leads 
to the formation of transient clusters and gels. This ef- 
fect increases with increasing attraction. The diffusion 
at short times is not influenced by collisions so that in 
the absence of attraction Dg is equal to the free diffu- 
sion of the particles. The decrease of Dg with increasing 
attraction is caused solely by cluster formation and can 
be calculated from the cluster size distribution. Since 
for any finite interaction there is a finite fraction of free 
particles and clusters, Ds only becomes zero when the 
interaction is infinitely strong. The decrease of Dg with 
increasing interaction strength is mainly determined by 
nin/F, see Fig. [9] and it is dominated by the decrease of 
the sol fraction for strong attraction when to„ is close to 
unity. 

The subsequent decrease of the diffusion coefficient 
from Dg to Di at long times is caused by crowding. An 
important observation is that the slowing down caused 
by crowding is independent of the attraction strength for 
a = 1. The reason is that by definition for reversible 
DLCA the reversibility is expressed as soon as collisions 
occur and the memory of the connectivity is lost. Conse- 
quently, diffusion is slowed down by attraction, but does 
not become zero as long as the attraction is finite. The 
situation is different for reversible RLCA in which case 
the collisions occur between long-lived clusters and be- 
tween the clusters and the percolating network. For given 
attraction strength Dg is the same for reversible DLCA 
and RLCA, but for RLCA the ratio Di/Dg decreases with 
decreasing bond formation probability. This situation 
will be explored elsewhere. 

The results presented here apparently contradict 
molecular dynamics simulations that showed critical 
slowing down at finite interaction strength f23, Is^ . 
The main difference between our simulation method and 
molecular dynamics is that rigid bonds are formed so that 
only cluster motion is possible. In the molecular dynam- 
ics simulations bound particles are still allowed to diffuse 
freely as long as no bonds are broken. The implication is 
that Dg is equal to the free particle diffusion independent 
of the interaction strength. Consequently Di is faster and 
Di is even larger than Dhg for weak attraction. We have 
included bond flexibility in the simulations with finite 
interaction strength by allowing free diffusion for bound 
particles as long as it does not lead to bond breaking. 
Details of these simulations will be reported elsewhere. 
Here, we only mention that including bond flexibility in- 
creases Di. It is clear, that allowing more freedom for 
movement does not lead to a critical slowing down of the 
diffusion at a finite interaction range so that the apparent 
contradiction with molecular dynamics simulations per- 
sists. But, a close look at the molecular dynamics simu- 
lation results shows that they can also be interpreted in 
terms of a power law decrease. For instance, Puertas et 
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FIG. 10: Long time self diffusion coefficient as a fuction of {(pp) 
for 2 different values of cj>p (0.4265 (circles) and 0.4519 (trian- 
gles)) from [H^]. The solid lines represents Di oc exp(— 50<^p). 



al. ^] simulated hard spheres with a short range deple- 
tion interaction caused by the addition of polymers. In 
Fig. [To] we have replotted Di as a function of the poly- 
mer volume fraction {(jjp oc —u) for two different particle 
volume fractions. The authors interpreted the data in 
terms of EqlJl At each concentration the smallest values 
of Di deviated from this expression, which was attributed 
to "hopping". It is clear, however, from Fig. [5] that the 
decrease of In Di at large (j>p can also be described as a 
power law even for the smallest value of Di implying that 
Di = only at infinite attraction. 

Very recently, Zaccarelli et al. [s^ made a detailed 
study of the dynamics in attractive hard sphere systems 
with limited valence (3 and 4) using molecular dynamics 
simulations. They observed that for strong attraction the 
variation of In Di could be described by a power law in 
terms of u at least for (j) < 0.55 implying that arrest only 
occurred at infinite attraction. On the other hand they 
could describe Di as a function of (p in terms of Eq. [T] (j)c 
was almost constant when increasing the attraction and 
only for strong attraction did they find a weak decrease 
of 0c, but the< extrapolation was uncertain in this case. 
The authors used the expression reversible gel for systems 
that arrested at u — > — oo (T — > 0) and glass (attractive 
or repulsive) for systems that arrested at (j) —^ (pc- The 
expression attractive glass was introduced to describe the 
arrest that occurs for a given volume fraction at a finite 
attraction energy. It is clear from the present study that 
attractive glasses in this sense are only formed for irre- 
versible aggregation. 
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V. CONCLUSION 

Reversible cluster aggregation of hard spheres leads 
to cqiiilibrium systems containing transient clusters and, 
above a critical interaction strength, a transient percolat- 
ing network. If the aggregation is diffusion limited and 
the bonds are rigid, then the effect of attraction on D; is 
decoupled from the effect of crowding. The latter is equal 
to that of non-interacting hard spheres, while the former 
is fully determined by the cluster size distribution. The 
self diffusion coefiicient of the spheres decreases with in- 
creasing attraction, but becomes zero only for irreversible 



aggregation, contrary to predictions from mode coupling 
theory. Therefore attractive glasses in the sense of sys- 
tems that are dynamically arrested by a finite interaction 
energy do not exist. 



Acknowledgments 

This work has been supported in part by a grant from 
the Marie Curie Program of the European Union num- 
bered MRTN-CT-2003-504712. 



[1] S. A. Safran, I. Webman, and G. S. Grest, Phys. Rev. A 

32, 506 (1985). 
[2] A. Rotenberg, J. Chem. Phys. 43, 1198 (1965). 
[3] B. J. Alder, D. A. Young, and M. A. Mark, J. Chem. 

Phys. 56, 3013 (1972). 
[4] D. Henderson, W. G. Madden, and D. D. Fitts, J. Chem. 

Phys. 64, 5026 (1976). 
[5] A. L. R. Bug, S. A. Safran, G. S. Grest, and I. Webman, 

Phys. Rev. Lett. 55, 1896 (1985). 
[6] S. C. Netemeyer and E. D. Glandt, J. Chem. Phys. 85, 

6054 (1986). 

[7] G. A. Chapela, S. E. Martinez-Casas, and C. Varea, J. 

Chem. Phys. 86, 5683 (1987). 
[8] Y. C. Chiew and Y. H. Wang, J. Chem. Phys. 89, 6385 

(1988). 

[9] E. Dickinson, C. Elvingson, and S. R. Euston, J. Chem. 
Soc., Faraday Trans. 2 85, 891 (1989). 
[10] D. M. Heyes, J. Chem. Soc. Faraday Trans. 87, 3373 
(1991). 

[11] L. Vega, E. de Miguel, L. F. RuU, G. Jackson, and I. A. 

McLure, J. Chem. Phys. 96, 2296 (1992). 
[12] A. Lomakin, N. Asherie, and G. B. Benedek, J. Chem. 

Phys. 104, 1646 (1996). 
[13] N. Asherie, A. Lomakin, and G. B. Benedek, Phys. Rev. 

Lett. 77, 4832 (1996). 
[14] N. V. BrilUantov, and J. P. Valleau, J. Chem. Phys. 108, 

1115 (1998). 

[15] J. R. ElUott and L. Hu, J. Chem. Phys. 110, 3043 (1999). 

[16] F. Del Rio, E. Avalos, R. Espindola, L. F. Ruh, G. Jack- 
son, and S. Lago, Mol. Phys. 100, 2531 (2002). 

[17] S. B. Kisclev, J. F. Ely, L. Lue, and J. R. EUiott, Fluid 
Phase Equihb. 200, 121 (2002). 

[18] G. Foffi, K. A. Dawson, S. V. Buldyrev, F. Sciortino, 

E. Zaccarelli, and P. Tartaglia, Phys. Rev. E 65, 050802 
(2002). 

[19] E. ZaccareUi, G. Foffi, K. A. Dawson, S. V. Buldyrev, 

F. Sciortino, and P. Tartaglia, Phys. Rev. E 66, 041402 
(2002). 

[20] S. Katsura, Phys. Rev. 115, 1417 (1959). 
[21] P. Orea, Y. Duda, and J. Alejandre, J. Chem. Phys. 118, 
5635 (2003). 

[22] P. Orea, Y. Duda, V. C. Weiss, W. Schrder, and J. Ale- 
jandre, J. Chem. Phys. 120, 11754 (2004). 

[23] G. Foffi, F. Sciortino, E. Zaccarelli, and P. Tartaglia, J. 
Phys.: Condens. Matter 16, S3791 (2004). 

[24] E. Zaccarelli, F. Sciortino, and P. Tartaglia, J. Phys.: 



Condens. Matter 16, S4849 (2004). 
[25] G. Foffi, E. Zaccarelli, S. V. Buldyrev, F. Sciortino, and 

P. TartagUa, J. Chem. Phys. 120, 8824 (2004). 
[26] G. Foffi, C. De Michelc, F. Sciortino, and P. Tartaglia, 

Piiys. Rev. Lett. 94, 078301 (2005). 
[27] G. Foffi, C. De Michele, F. Sciortino, and P. Tartaglia, J. 

Chem. Phys. 122, 224903 (2005). 
[28] D. L. Pagan and J. D. Gunton, J. Chem. Phys. 122, 

184515 (2005). 

[29] S. Babu, J.-C. Gimel, and T. Nicolai, J. Chem. Phys. 

125, 184512 (2006). 
[30] M. A. Miller and D. Frenkel, Phys. Rev. Lett. 90, 135702 

(2003) . 

[31] M. A. Miller and D. Frenkel, J. Chem. Phys. 121, 535 

(2004) . 

[32] S. Diez Orritc, S. StoU, and P. Schurtenberger, Soft Mat- 
ter 1, 364 (2005). 

[33] P. Mcakin, J. Chem. Phys. 83, 3645 (1985). 

[34] M. Kolb, J. Phys. A 19, L263 (1986). 

[35] W. Y. Shih, I. A. Aksay, and R. Kikuchi, Phys. Rev. A 
36, 5015 (1987). 

[36] W. Y. Shih, J. Liu, W.-H. Shih, and I. A. Aksay, J. Stat. 
Phys. 62, 961 (1991). 

[37] M. D. Haw, M. Sievwright, W. C. K. Poon, and P. N. 
Puscy, Adv. Colloid Interface Sci. 62, 1 (1995). 

[38] J.-M. Jin, K. Parbhakar, L. H. Dao, and K. H. Lee, Phys. 
Rev. E 54, 997 (1996). 

[39] T. Terao and T. Nakayama, Phys. Rev. E 58, 3490 
(1998). 

[40] T. Terao and T. Nakayama, Prog. Theor. Phys. Suppl. 

138, 354 (2000). 
[41] J. C. Gimel, T. Nicolai, and D. Durand, Eur. Phys. J. E 

5, 415 (2001). 

[42] G. Odriozola, A. Schmitt, A. Moncho-Jorda, J. Callejas- 

Fernandcz, R. Martincz-Garcia, R. Leone, and R. 

Hidalgo- Alvarez, Phys. Rev. E 65, 031405 (2002). 
[43] J. C. Gimel, T. Nicolai, and D. Durand, Phys. Rev. E 

66, 061405 (2002). 
[44] G. Odriozola, A. Schmitt, J. Callejas-Fernandez, R. 

Martinez-Garcia, R. Leone, and R. Hidalgo-Alvarez, J. 

Phys. Chem. B, 107, 2180 (2003). 
[45] S. Babu, M. Rottereau, T. Nicolai, J. C. Gimel, and D. 

Durand, Eur. Phys. J. E 19, 203 (2006). 
[46] R. J. Speedy, Mol. Phys. 62, 509 (1987). 
[47] H. Sigurgeirsson and D. M. Heyes, Mol. Phys. 101, 469 

(2003). 



9 



[48] I. Moriguchi, J. Chem. Phys. 106, 8624 (1997). 
[49] B. Cichocki and K. Hinsen, Physica A 187, 133 (1992). 
[50] B. Doliva and A. Heuer, Phys. Rev. E 61, 6898 (2000). 
[51] F. Sciortino and P. Tartaglia, Adv. Phys. 54, 471 (2005). 
[52] W. van Megen, T. C. Mortensen, S. R. Wilhams, and J. 

Miiller, Phys. Rev. E 58, 6073 (1998). 
[53] T. Voigtmann, A. M. Puertas, and M. Fuchs, Phys. Rev. 

E 70, 061506 (2004). 
[54] A. M. Puertas, M. Fuchs, and M. E. Gates, Phys. Rev. 

E 67, 031406 (2003). 

[55] S.Babu, J. C. Gimel and T. Nicolai a rXiv:0705.1266 ^yl 

[cond- mat. soft]. 
[56] K. A. Dawson, G. Fofii, M. Fuchs, W. Gotze, F. Sciortino, 

M. Sperl, P. Tartaglia, T. Voigtmann and E. Zaccarelh, 

Phys. Rev. E 63 011401 (2000). 
[57] K. N. Pham, S.U. Egelhaaf, P. N. Pusey, and W. G. K. 

Poon, Phys. Rev. E 69, 011503 (2004). 
[58] M. Rottereau, J. G. Gimel, T. Nicolai, and D. Durand, 

Eur. Phys. J. E 15, 133 (2004). 



[59] E. Zaccarelh, I. Saika-Voivod, S. V. Buldyrev, A. J. 

Moreno, P. Tartaglia, and F. Sciortino, J. Ghem. Phys. 

124, 124908 (2006). 
[60] E. Zaccarelh, S. V. Buldyrev, E. La Nave, A. J. Moreno, I. 

Saika-Voivod, F. Sciortino, and P. Tartaglia, Phys. Rev. 

Lett. 94, 218301 (2005). 
[61] E. Del Gado and W. Kob, Europhys. Lett. 72, 1032 

(2005). 

[62] F. Sciortino, P. Tartaglia, and E. ZaccareUi, J. Phys. 

Chem. B (2005). 
[63] A. Gonigho, L. de Arcangelis, E. del Gado, A. Fierro, and 

N. Sator, J. Phys. Gondens. Matter 16, S4831 (2004). 
[64] M. G. Noro and D. Frenkel, J. Ghem. Phys. 113, 2941 

(2000). 

[65] T. Kihara, Rev. Mod. Phys. 25, 831 (1953). 
[66] A. M. Puertas, M. Fuchs, and M. E. Gates, J. Phys. 
Ghem. B 109, 6666 (2005). 



